Transverse-Mass Effective Temperature in Heavy-Ion Collisions from AGS to SPS 
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Transverse-mass spectra in Au+Au and Pb+Pb collisions in incident energy range from 2A to 
160A GeV are analyzed within the model of 3-fluid dynamics (3FD). It is shown that dynamical 
description of freeze-out, accepted in this model, naturally explains the incident energy behavior of 
inverse-slope parameters of these spectra observed in experiment. Simultaneous reproduction of the 
inverse-slopes of all considered particles (p, n and K) suggests that these particles belong to the 
same hydrodynamic flow at the instant of their freeze-out. 
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Experimental data on transverse-mass spectra of kaons 
produced in central Au+Au l] or Pb+Pb Q collisions 
reveal peculiar dependence on the incident energy. The 
inverse-slope parameter (so called effective temperature 
T) of these spectra at mid rapidity increases with incident 
energy in the energy domain of BNL Alternating Gradi- 
ent Synchrotron (AGS) and then saturates at energies of 
CERN Super Proton Synchrotron (SPS). In Refs. @, 
it was assumed that this saturation is associated with the 
deconfinement phase transition. This assumption was in- 
directly confirmed by the fact that microscopic transport 
models, based on hadronic degrees of freedom, failed to 
reproduce the observed behavior of the kaon inverse slope 
[SHI- Hydrodynamic simulations of Ref. Q succeeded to 
describe this behavior. However, in order to reproduce it 
these hydrodynamic simulations required incident-energy 
dependence of the freeze-out temperature which almost 
repeated the shape of the corresponding kaon effective 
temperature. This happened even in spite of using equa- 
tion of state (EoS) involving the phase transition into 
quark-gluon plasma (QGP). This way, the puzzle of kaon 
effective temperatures was just translated into a puzzle of 
freeze-out temperatures. Moreover, results of Ref. im- 
ply that peculiar incident-energy dependence of the kaon 
effective temperature may be associated with dynamics 
of freeze-out. 

In this paper we would like to present calculations of 
effective temperatures within the 3FD model 0, H, |T(| 
which is suitable for simulating heavy-ion collisions in the 
range from AGS to SPS energies. We perform our sim- 
ulations [1, [l(| with a simple, hadronic EoS [HJ which 
involves only a density dependent mean field providing 
saturation of cold nuclear matter at normal nuclear den- 
sity and with the proper binding energy. The 3FD model 
with the intermediate EoS turned out to be able to rea- 
sonably reproduce a great body of experimental data [|[ 
in a wide energy range from AGS to SPS. In particu- 
lar, transverse-mass spectra of protons were reproduced. 
This was done with the unique set of model parameters 
summarized in Ref. 0]. Problems were met only in de- 
scription of transverse flow The directed flow re- 
quires a softer EoS at top AGS and SPS energies (in 
particular, this desired softening may signal occurrence 



of the phase transition into QGP). Similar softening is 
needed for reproduction of recent data on rapidity distri- 
butions of net-baryon number in central Pb+Pb collision 
at energies 20A-80A GeV [nj]. 

The transverse-mass spectra are most sensitive to the 
freeze-out parameters of the model. In fact, inverse 
slopes of these spectra represent a combined effect of 
the temperature and collective transverse flow of expan- 
sion. Fig. [T] demonstrates these important interplay. 
Had it been only the effect of thermal excitation, inverse 
slopes for different hadronic species would approximately 
equal. The collective transverse flow makes them differ- 
ent. These two effects partially compensate each other: 
the later freeze-out occurs, the lower temperature and 
the stronger collective flow arc. Nevertheless, transverse- 
mass spectra turn out to be sensitive to the instant of 
the freeze-out. 

3FD results for inverse-slope parameters of transverse- 
mass spectra of kaons, pions and protons produced in 
central Au+Au and Pb+Pb collisions are presented in 
Fig. [TJ The inverse slopes T were deduced by fitting the 
calculated spectra by the formula 
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where hit and y are the transverse mass and rapidity, re- 
spectively. Though the purely exponential fit with A = 
does not always provide the best fit of the spectra, it al- 
lows a systematic way of comparing spectra at different 
incident energies. In order to comply with experimental 
fits at AGS energies (and hence with displayed experi- 
mental points), we also present results with A = — 1 for 
pions and with A = 1 for protons. These results are ob- 
tained with precisely the same EoS and set of parameters 
(friction, freeze-out and formation time) as those used in 
Ref. [||, which was found to be the best for other ob- 
servables. No special tuning was done to reproduce these 
effective temperatures. 

Numerical problems, discussed in Ref. Q, prevented 
us from simulations at RHIC energies. Already for the 
central Pb+Pb collision at the top SPS energy the code 
requires 7.5 GB of (RAM) memory. At the top RHIC en- 
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ergy, required memory is three order of magnitude higher, 
which is unavailable in modern computers. 
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FIG. 1: (Color online) Inverse-slope parameters of transverse- 
mass spectra of kaons, pions and protons at mid rapidity pro- 
duced in central Au+Au and Pb+Pb collisions as a function 
of invariant incident energy. Solid lines correspond to purely 
exponential fit (A = 0, see Eq. (HJ), while dashed lines present 
results with A = — 1 for pions and with A = 1 for protons. Ex- 
perimental data are from Refs. HHQ. 

As seen from Fig. [TJ reproduction of effective tempera- 
tures is quite reasonable. Moreover, the pion and proton 
effective temperatures also reveal saturation at SPS en- 
ergies, if they are deduced from the purely exponential fit 
with A = 0. It is important that it is achieved with a sin- 
gle freeze-out parameter £f rz = 0.4 GeV/fm 3 , the critical 
freeze-out energy density, which is the same for all con- 
sidered incident energies above 2A GeV, both for chemi- 
cal and thermal freeze-out. Only for smaller energies we 
used smaller values: Ef rz (2A GeV) = 0.3 GeV/fm 3 and 
e irz (lA GeV) = 0.2 GeV/fm 3 . In order to clarify why 
this happens, let us turn to the 3FD freeze-out proce- 
dure, which is analyzed in Rcf. (l~5j in more detail. 

The freeze-out criterion we use is 



£ < £ fr; 



(2) 



where e = u^T^u^ is the total energy density of all 
three fluids in the proper reference frame, where the 
composed matter is at rest. This total energy density 
is defined in terms of the total energy-momentum ten- 
sor = Tg v + Tl 1 " + Tf v being the sum of energy- 
momentum tensors T£" of separate fluids (projectile-like, 
target-like and fireball ones) and the total collective 4- 
velocity of the matter = u v T^ v /u\T Xk u k . Note the 
latter definition is, in fact, an equation determining v} 1 . 
A very important feature of our freeze-out procedure is 



an anti-bubble prescription. The matter is allowed to be 
frozen out only if 

(a) either the matter is located near the boarder with 
vacuum (this piece of matter gets locally frozen out) 

(b) or the maximal value of the total energy density in 
the system is less than £f rz 



maxe < £f rz 



(3) 



(the whole system gets instantly frozen out). 

In the 3FD model this freeze-out simultaneously ter- 
minates both chemical and kinetic processes. 

Before the instant of the global freeze-out, cf. j3|, 
the freeze-out remove matter from the surface of the hy- 
drodynamically expanding system. This removed matter 
gives rise to observable spectra of hadrons. This kind of 
freeze-out is similar to the model of "continuous emis- 
sion" proposed in Rcf. There the particle emission 
occurs from a surface layer of the mcan-free-path width. 
In our case the physical pattern is the same, only the 
mean free path is shrunk to zero. 
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FIG. 2: (Color online) Actual average freeze-out energy den- 
sity in central (zero impact parameter) Pb+Pb collisions as a 
function of invariant incident energy. 

Condition @ ensures only that the actual freeze-out 
energy density (let us call it £ ou t), at which the freeze-out 
actually occurs, is less than £f rz . Therefore, £f rz can be 
called a "trigger" value of the freeze-out energy density. 
As explained in Ref. a natural value of this actual 
freeze-out energy density is £ ou t ~ £s/2, i.e. at that 
the middle of the fall from the near-surface value of the 
energy density, s s , to zero. To find out the actual value of 
£out; w e have to analyze results of a particular simulation. 
In our previous paper Q we have performed only a rough 
analysis of this kind. This is why in the main text of 
Ref. [1[ we mentioned the value of approximately 0.2 
GeV/fm 3 for £ out and in appendix explained how the 
freeze-out actually proceeded. [In terms of Ref. @) ( e frz[il 
and crti) our present quantities are £f rz = Ejwi] an< ^ 
£out = £frz[i]-] Results of more comprehensive analysis 
for central (b = 0) Pb+Pb collisions are presented in Fig. 
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which shows the £ ou t value averaged over space-time 
evolution of the collision: (e out )- As seen, (e ut) reveals 
saturation at the SPS energies, very similar to that in 
effective temperatures in Fig. [T] This happens in spite 
of the fact that our freeze-out condition involves only a 
single constant parameter £f rz . 

The "step-like" behavior of (e ou t) is a consequence of 
the freeze-out dynamics, as it was demonstrated in Ref. 
[l5l ]. At low (AGS) incident energies, the energy density 
achieved at the border with vacuum, e s , is lower than 
£f rz . Therefore, the surface freeze-out starts at lower en- 
ergy densities. It further proceeds at lower densities up to 
the global freeze-out because the freeze-out front moves 
not faster than with the speed of sound, like any pertur- 
bation in the hydrodynamics. Hence it cannot overcome 
the supersonic barrier and reach dense regions inside ex- 
panding system. With the incident energy rise the energy 
density achieved at the border with vacuum gradually 
reaches the value of £f rz and then even overshoot it. If 
the overshoot happens, the system first expands with- 
out freeze-out. The freeze-out starts only when e s drops 
to the value of £f rz . Then the surface freeze-out occurs 
really at the value e s ~ £f rz and thus the actual freeze- 
out energy density saturates at the value (e ou t) ~ £frz/2. 
This freeze-out dynamics is quite stable with respect to 
numerics 11511. 
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FIG. 3: (Color online) Average temperature (Tf rz , see scale 
at the l.h.s. axis) , transverse velocity (j3 = vt / c, see scale at 
the r.h.s. axis) and baryon density over the normal nuclear 
density (rif rz /no, see scale at the r.h.s. axis) at the freeze-out 
in central Au+Au (at AGS energies, b = 2 fm) and Pb+Pb (at 
SPS energies, b — 2.5 fm) collisions as a function of invariant 
incident energy. Arrows point to the proper axis scale. 

Fig. [3] presents average temperatures, transverse ve- 
locities and baryon densities achieved at the freeze-out in 
central collisions. The freeze-out temperature Tf rz has a 
similar "step-like" behavior. However, its absolute values 
are essentially lower than effective temperatures in Fig. 
[TJ This fact once again illustrates that inverse slopes rep- 



resent a combined effect of the temperature and collective 
transverse flow of expansion associated with the collec- 
tive transverse velocity /3. Note that (3 even slightly de- 
creases with incident energy. At SPS energies the freeze- 
out temperatures in Fig. [3] are noticeably lower than 
those deduced from hadron multiplicities in the statis- 
tical model [13, [H|]. The reason for this is as follows. 
Whereas the statistical model assumes a single uniform 
fireball, in the 3FD simulations at the late stage of the 
evolution the system effectively consists of several "fire- 
balls" : two (one baryon-rich and one baryon-free) fire- 
balls at lower SPS energies and three (two baryon-rich 
and one baryon-free) fireballs at top SPS energies [l5l |. 
Therefore, whereas high multiplicities of mesons and an- 
tibaryons are achieved by means of high temperatures 
in the statistical model, the 3FD model explains them 
by an additional contribution of the baryon-free fireball 
at a lower temperature. In particular, this is the reason 
why two different freeze-out points (chemical and kinetic 
ones) arc not needed in the 3FD model. The freeze-out 
baryon density nf rz exhibits a maximum at incident ener- 
gies of £i a b = 10A-30^4 GeV which are well within range 
of the planned FAIR in GSI. This observation agrees with 
that deduced from the statistical model [lj| , even baryon 
density values in the maximum arc similar to those pre- 
sented in Ref. fl9l |. 

Returning to the question if the considered " step- like" 
behavior of effective temperatures is a signal of phase 
transition into QGP, wc should admit that this is not 
quite clear as yet. It depends on the nature of the freeze- 
out parameter £f rz = 0.4 GeV/fm 3 which should be fur- 
ther clarified. EoS is not of prime importance for this 
behavior. The only constrain on the EoS is that it should 
be in some way reasonable. Moreover, our preliminary 
results indicate that a completely different EoS [2(| with 
1st order phase transition to QGP still reasonably re- 
produces this " step- like" behavior even in spite of that it 
fails to describe a large body of other data. This happens 
because the same freeze-out pattern is accepted there. 

In fact, EoS is just the pressure as a function of baryon 
and energy densities: P(ub, £)■ In this calculation we 
used hadronic EoS [1, However, this is just an 

interpretation of the function P(jib,£), which we use, 
in hadronic terms. Moreover, our EoS is too soft at 
high densities to be matched with even heavy-quark bag- 
model EoS [2(| or quasiparticle fits to lattice QCD data 
[2l| in order to construct the 1st order phase transi- 
tion. [The hadronic and quark pressures as functions 
of the baryon chemical potential should have a crossing 
point in order to construct the 1st order phase transi- 
tion.] Therefore, it would not be surprising if the same 
EoS can be reinterpreted also in terms of very smooth 
cross-over phase transition to quark-gluon matter. 

This hydrodynamic explanation of the considered 
"step-like" behavior of effective temperatures together 
with failure of kinetic approaches implies that a heavy nu- 
clear system really reveals a hydrodynamic motion dur- 
ing its expansion. Simultaneous reproduction of inverse- 
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slopes of all considered particles (p, tt and K) implies that 
these particles belong to the same hydrodynamic flow at 
the instant of their freeze-out. An indirect support of 
this conjecture is the recent success of the GiBUU model 
[22| in reproduction of kaon inverse-slopes. That was 
achieved by taking into account three-body interactions, 
which essentially increased the equilibration rate. 
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